Scaling in soft spheres: fragility invariance on the repulsive potential softness 
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We address the question of the dependence of the fragility of glass forming supercooled liquids on 
the "softness" of an interacting potential by performing numerical simulation of a binary mixture 
of soft spheres with different power n of the interparticle repulsive potential. We show that the 
temperature dependence of the diffusion coefficients for various n collapses onto a universal curve, 
supporting the unexpected view that fragility is not related to the hard core repulsion. We also find 
that the configurational entropy correlates with the slowing down of the dynamics for all studied n. 

PACS numbers: 61.20.Ja, 02.70.Ns, 64.70.Pf 
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When a liquid is cooled below its melting temperature, 
if crystallization does not take place, it becomes super- 
cooled. In this supercooled region, the viscosity increases 
by more than 15 order of magnitude in a a small T-range. 
When the viscosity rj reaches a value of about 10 13 Poise 
the liquid can be treated as an amorphous solid, i.e., a 
glass Ql [3] and the corresponding temperature is 
defined as the glass transition temperature (labelled T g ). 

The T-dependence of the viscosity r\ differs for different 
glass formers. Angell has proposed a classification based 
on the behaviour of r/(T). Glasses are said to be frag- 
ile if they show large deviations from an Arrhenius law 
(i](T) oc exp[2?/T]) or strong otherwise The fragility 
m of a glass forming liquid can be quantified by the slope 
of log tj(T) vs Tg/T, evaluated at T g , i.e. as 



dlogry 
d(T g /T) 



\T=Tg 



(1) 



While the original definition of fragility is based on a 
purely dynamic quantity, correlation between m and 
other physical properties of glass forming liquids, both 
with dynamic and thermodynamic properties, have been 
reported . Recently, a remarkable correlation with vibra- 
tional properties of the glass state has been discovered 
. One of the main challenges in the physics of super- 
cooled liquid and glasses is to understand the connection 
between dynamical properties of the liquid close to the 
glass transition, i.e. the fragility, and microscopic prop- 
erties. Is the fragility most affected by the steepness of 
the repulsive potential or by the inter particle attraction? 
Is it controlled by other properties of the interaction po- 
tential? In the present letter we address this question 
calculating numerically the fragility of several models for 
liquids, differing only in the softness of the repulsive po- 
tential. We aim at understanding whether changing the 
softness of the repulsive potential the fragility changes 
accordingly. We show more generally that the diffusion 
coefficient D can be scaled on a universal master curve 
by changing the softness of the repulsive potential. This 



implies that the fragility does not depend on the softness 
of the interaction potential. We complement this dynam- 
ical study with the evaluation of the configurational en- 
tropy to check the validity of the Adam-Gibbs fll lil ioL 
relation. 

In this Letter we consider a simple glass former, a 80 : 
20 binary mixture of N = 1000 soft spheres [Til [la. fl3| . 
that is an ensemble of spheres interacting via the follow- 
ing potential 



V a p{r)=U aP (?f) n 



(2) 



where a, (3 S A,B, gaa — 1-0, oab — 0.8, obb = 0.88, 
£,4,4 = 1-0, £ab — F5, €bb = 0.5 and n is a parameter 
by which is possible to tune the "softness" of the inter- 
action [lJLThis interaction potential is a Kob- Andersen 
potential |l5j in which the attractive part of the poten- 
tial has been dropped. In particular we investigate the 
values n = 6, 8, 12, 18. 

This choice for the binary mixture is motivated by the 
fact that such a system is not prone to crystallization, 
that is it can be easily supercooled below its melting 
temperature . Still, for n < 6, crystallization takes place 
within the simulation time, determining a lower limit to 
the range of investigated n values. Reduced units will be 
used in the following, length will be in units of uaa, en- 
ergy in units of caa and time in units of (Ma AA / '^aa) 1 ^ 
where M is the mass of all particles. In physical units, 
assuming the atom A is Argon, the units are a length of 
3.4A, an energy of \2QKks and a time of 2.1bps. 

The self-similar nature of the soft-sphere potential cou- 
ples T and V . It can be shown that all thermodynamic 
properties depend on the quantity TVn [la. Dynamic 
properties can also be scaled accordingly |ljj. Hence, 
it is sufficient to quantify the T-dependence or the In- 
dependence of any observable to fully characterize the 
behavior of the system. As a consequence the fragility 
does not change upon changing the density of the soft 
binary mixture. 
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FIG. 1: T-dependence of the diffusio n coefficient for all n 
investigated. D is measured in units y <j\ a £aa/M ■ 



Figure ^ shows the T-dependence of the diffusion co- 
efficients, evaluated from the long time limit behavior of 
the mean square displacement, for all n investigated and 
covering a window of about four order of magnitudes. 

In the attempt of compare the n-dependence of the 
diffusion coefficient, we report in Fig. [21 the data as a 
function of T n /T, where T n is chosen in such a way to 
maximize the overlap between data of different n, i.e. to 
collapse all data onto a single master curve. Figure 
shows that all curves can be successfully scaled onto the 
master curve T> choosing a proper set of scaling parame- 
ters T n (whose n-dependence is plotted in the inset of this 
Figure). The very good quality of the resulting master 
curve 



D(T) = V(T/T n ). 



(3) 



suggests that the n-dependence enters only via a rescaling 
of the temperature The remarkable consequence 

of latter result is that the fragility of the system does not 
depend on the repulsive interaction potential. In fact 
according to Eq. |J3J| and from the definition of liquid's 
fragility m given in Eq. ijTJl. assuming D oc t _1 we get: 



f 



T g (n) 

T n V[T g {n)/T n 



dV(x) 
dx 



U=T fl (n)/T„ 



(4) 



where T g (n) is the glass transition temperature for the 
system with softness n, which can be defined as the tem- 
perature at which diffusivity reaches an arbitrary small 
value lO* E3, i.e. 



log D[T g (n)} = -logX> 



T g (n) 



= K 



(5) 



Eq0] shows that the fragility index m is a function 



only of the scaled variable 



T g (n) 



and hence, as far as the 



FIG. 2: Master curve for diffusion coefficient D calculated 
rescaling the temperature by an n-dependent quantity. Full 
line is a Vogel-Fulcher-Tamman fit to all points lying in the 
landscape- influenced region |20l l2ll| . which for the present 
model corresponds to £ < 0.78 with £ = T/T», , i.e. £>(£) = 
Aexp[B / — fo)]. The fitting parameters In A, B and £o are 
respectively —2.46,-0.857,0.209. From this master curve one 
can also estimate the dynamic fragility, that is about 130.4. 
Dashed line is a fit to all points lying in the interval 0.328 < 
£ < 0.833, to a power law [^lEJEl], i.e. £>(£) = A m (£-£c) 7 , 
where A m = -2.41,£ c = 0.309,7 = 1-88. The inset shows the 
n-dependence of the scaling parameter T n . 



scaling reported in Fig. |2 keeps holding even at tem- 
peratures lower than the one we are able to equilibrate, 
the dynamic fragility m is independent of n as well. By 
fitting the master curve to a Vogel-Tamman-Fulcher fit, 
as shown in Fig. an estimate of T j! = 10*" can be 
calculated, resulting into a estimation of m « 130. This 
figure should be compared with the value m = 81 for 
o-terphenyl (OTP), that is a typical fragile liquid and 
rn = 20 for the prototypical strong glass the liquid silica 
(Si0 2 ). 

For completeness, we report also in Fig. [21 a fit of 
the master curve according to the prediction of mode- 
coupling theory, which has been shown to be consistent 
with numerical data for several models in the weak su- 
percooling region. A best fit procedure requires the ex- 
clusion of the low T points, for which deviations from the 
power-law fit are observed |2(| . 

Recently, evidence has been presented that ki- 
netic fragility strongly correlates with thermodynamic 
fragility [2j . In this respect, it is worth looking if the scal- 
ing observed in dynamical properties has a counterpart 
in thermodynamic properties. In particular, we evaluate 
the configurational entropy for the system, within the 
potential energy landscape framework as discussed in de- 
tails in Refs.[2imm|3jm|33!]. In brief, we estimate 
S c as difference between the liquid entropy (calculated 
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FIG. 3: Test of the Adam-Gibbs relation based on the con- 
figurational entropy values calculated subtracting the anhar- 
monic and harmonic entropies to the total entropy (see |2£l ] 
for details). The fitting parameters In Aag and Bag are: for 
n = 6, -1.112, -1.584; for n = 8, -1.303, -2.105 ; for n = 12, 
-1.271, -3.481; for n = 18, -1.657, -4.002. TS C is measured 
in units of €aa- 

via thermodynamic integration from the ideal gas) and. 
of the vibrational entropy (calculated via thermodynamic 
integration, including anharmonic corrections, from the 
very low temperature harmonic dynamics of the disor- 
dered solid associated to the liquid configuration). We 
then focus on the ability of the Adam-Gibbs (AG) rela- 
tion — which states that 

D(T)=A AG e%£, (6) 

— of modelling the temperature dependence of D. Fig. [3] 
shows the AG plot for the studied n values. For all n, a 
satisfactory linear representation of log(D) vs. 1/TS C (T) 
is observed. As discussed in more details in Ref . . the 
simultaneous validity of the VTF description of D and of 
the AG relation requires the identity of the kinetic and 
thermodynamic fragilities. In this respect, the indepen- 
dence of n discussed above for the case of kinetic fragility 
carries on also to thermodynamic fragility. 

A remarkable consequence of the validity of the AG 
relation (Eq. EJ), associated to the scaling with n of D 
(Eq|3J) is that the configurational entropy can be written 
as 

S C (T) = S {n)T(T/T n ) (7) 

where F(x) is a scaling function and Sq(ji) = Bag/T u . 
To support such proposition, we show in Figure 0] S c 
multiplied by the factor Bag/Tu as a function of T/T n , 
were T n are the values for which D scaling is recovered 
(inset of FigJ^Jl. Again, the quality of the data collapse 
stresses the validity of the scaling with n. 



FIG. 4: Master plot of the configurational entropy (see 
Eq. J7J). The inset shows the temperature dependence of 
the configurational entropy. 



To conclude the relevant result that has been shown in 
this Letter is that in the case of soft, sphere potentials, 
the dynamic fragility is independent on the power n of 
the short range repulsion. This conclusion is based on 
the hypothesis that the scaling observed in the range of 
T where simulations are feasible extends also to lower 
temperatures, down to the glass transition temperature. 
Indeed, a particular effort has been made to equilibrate 
configurations to temperatures lower than the MCT tem- 
perature, where dynamical processes different from the 
ones captured by MCT are active. If the scaling is in- 
deed valid, the results presented in this Letter strongly 
support the possibility, that contrary to our common un- 
derstanding, fragility in liquids is mostly controlled by 
other properties of the potential, more than by the hard 
core repulsion. Finally we note that one could be tempted 
to associate the fact that the diffusion coefficient data 
can be rescaled only by change the energy scale by T n 
to a simple overall rescaling of the landscape potential. 
The data in Fig^ suggest that this is not the case since 
S C (E) is not just scaling function of T/T n but it needs to 
be rescaled by a factor So(n) and hence the number of 
distinct basins explored at the same T/T n changes with 
n. A non-trivial compensation mechanism between the 
scaling of the static properties (S c ) and the scaling of the 
kinetic coefficient BAG{ n ) (defined in Eq. JBJ) on n must 
be present. 
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